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摘 要 : 小 流域 是 研究 小 微 尺度 水 文 水 资源 系统 演变 规律 的 理想 对 象 ,是 用 于 计算 河流 产 水 产 沙 
的 最 小 单元 ,是 水 文 及 水 土 流失 研究 与 管理 的 最 佳 地 域 久 度 。 通 过 遥感 技术 ,气候 模式 获得 降水 
数据 ,并 驱动 分 布 式 水 文 模型 ,模拟 和 预测 水 文 过 程 ,是 流域 水 文 水 资源 研究 的 必然 趋势 。 以 NO- 
AA-CPC-US 降水 作为 参照 ,在 美国 不 同 地 区 的 9 个 小 流域 ,评估 卫星 降水 产品 PERSIANN, PERSI- 
ANN-CDR,TRMM-3B42V7,GPM-IMERG , 雷达 降水 StageIV 以 及 气候 模式 ERA5 降水 产品 的 精度 ， 
并 用 这 7 种 降水 产品 驱动 CREST 分 布 式 水 文 模型 ,评估 了 7 种 降水 产品 的 水 文 模拟 效用 。 研 究 表 
明 : 各 降水 产品 与 NOAA-CPC-US 降水 吻合 程度 从 高 到 低 , 依 次 是 StagelV 雷达 降水 ,PERSIANN- 
CDR 和 GPM-IMERG 次 之 ,再 次 是 PERSIANN 和 ERA5, 最 后 是 TRMM-3B42V7。 各 降水 产品 在 美 
国 北 部 高 纬度 地 区 和 西部 山地 等 区 域 的 小 流域 降水 估算 精度 略 低 , 在 美国 中 部 ,南部 ,东部 的 小 流 
域 有 较 好 的 降水 精度 。 在 水 文 模拟 效用 评估 中 , 设 定 相同 率 定期 ,分 别 使 用 7 种 降水 产品 率 定 
CREST 模 型 参数 ,得 到 率 定 参数 集 后 ,在 相同 验证 期 对 流域 日 径流 过 程 进行 模拟 。 结 果 表 明 :NO- 
AA-CPC-US 和 Stage IV 雷达 降水 在 各 小 流域 水 文 模拟 中 效果 较 好 ,在 美国 北部 和 西部 地 区 ,使 用 
PERSIANN, PERSIANN- CDR, GPM- IMERG, ERAS 降水 进行 水 文 模拟 时 需要 谨慎 。TRMM- 


3B42V7 的 小 流域 水 文 模拟 效果 不 理想 。 
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降水 是 水 文 循环 的 主要 环节 之 一 ,也 是 开展 水 
文 模型 研究 的 主要 数据 输入 。 降 水 在 气候 ,气象 ， 
水 文 预报 ,水文 灾 害 等 领域 的 研究 ,监测 和 预测 中 
发 挥 着 重要 作用 ”。 受 气候 因素 和 下 垫 面 因 素 影 


面 雨 量 站 观测 手段 直接 ,精度 较 高 ,但 受 站 点 密度 
及 其 空间 分 布 影响 ,难以 给 出 复杂 多 变 的 降水 空间 
分 布 ;地 面 雷 达 测 雨 可 以 直接 探测 降水 的 空间 分 
布 ,具有 实时 跟踪 降水 中 心 和 空间 变化 的 能 力 ,但 


响 , 降 雨 的 形成 和 分 配 在 时 空 分 布 不 均匀 ,进而 又 
影响 了 陆地 径流 、 蒸 散发 .土壤 水 分 .地 下 水 等 流域 
水 循环 过 程 。 因 此 ,对 降水 的 准确 估计 和 大 范围 快 
速 获取 高 时 空 分 辩 紊 的 连续 降水 观测 数据 ,不仅 能 
够 促进 我 们 对 水 循环 的 理解 ,而 且 对 流域 水 文 过程 
模拟 以 及 预报 也 具有 重要 意义 。 目 前 ,获得 降水 数据 
的 方法 有 4 种 :地 面 雨量 站 观测 ,地 面 雷达 探测 降水 ， 


易 受 地 形迹 挡 , 雷 达 射 线 抬升 等 影响 *。 卫 星 膛 感 
测 雨 虽然 精度 不 如 地 面 观测 降水 精度 高 ,但 它 能 够 
大 范围 ,快速 .便捷 地 获取 具有 一 定时 空 精度 的 连 
续 降 水 数据 ,一 定 程度 上 能 弥补 地 面 观测 降水 的 不 
足 , 对 无 资料 和 少 地 面 观测 资料 地 区 是 一 种 有 效 的 
降水 数据 来 源 。 气 候 系统 模式 可 以 模拟 全 球 降 水 ， 
然而 受气 修复 杂 性 影响 ,其 模拟 结果 在 区 域 尺度 仍 


卫星 遥感 降水 和 气候 系统 模式 数值 模拟 降水 ,地 ”存在 较 大 不 确定 性 。 一 些 气候 模式 的 降水 模拟 
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数据 还 不 能 直接 用 于 水 文 效用 评估。 

近年 来 , 随 着 上 述 4 种 降水 数据 获取 技术 的 快 
速 发 展 , 以 及 机 器 学 习 算法 在 反 演 降水 中 的 应 用 ， 
涌现 出 大 量 不 同时 空 分 辩 率 的 降水 产品 “局 。 同 
时 ,诸多 学 者 针对 不 同 降水 产品 ,围绕 降水 精度 评 
佑 ,不同 时空 尺度 降水 变化 ,模拟 水 文 效用 等 方面 
开展 了 大 量 研究 ”7 ,基于 降水 产品 开发 的 全 球 或 大 
中 型 区 域 尺度 的 径流 模拟 和 实时 预报 系统 ,在 多 次 水 
文 灾害 事件 的 预报 和 再 现 中 得 到 了 积极 的 评价 “20。 

然而 ,前 人 研究 结论 已 发 现 降水 产品 目前 尚 存 
在 明显 的 空间 变异 性 ,卫星 ,雷达 等 降水 产品 在 不 
同 维度 地 区 的 精度 有 差别 ,而且 不 同 精度 的 降水 产 
品 又 会 导致 径流 模拟 精度 的 差别 , 即 多 源 降 水 产品 
在 不 同时 空 尺 度 的 水 文 应 用 中 仍 存在 明显 的 不 确 
定性 下。 小 流域 是 研究 小 微 尺 度 水 文 水 资源 系统 
演变 规律 的 理想 对 象 ,是 水 土 流失 治理 ,生态 系统 
修复 等 研究 工作 的 重要 地 理 单元 ,也 是 用 于 计算 河 
流产 水 产 沙 的 最 小 单元 ,是 水 文 及 水 土 流 失 研究 与 
管理 的 最 佳 地 域 尺度 。 近 年 来 ,学 者 在 小 流域 尺度 
基于 地 面 观 测 降水 围绕 产 流 系数 ,地表 径 流 模拟 ， 
土壤 水 运动 做 了 诸多 研究 工作 喇 汪 ,但 就 不 同时 空 
分 辩 率 的 降水 产品 在 小 流域 尺度 水 文 模 拟 和 预报 
中 的 能 力 还 有 待 深入 研究 。 

鉴于 此 ,本 文选 取 分 布 在 美国 9 个 不 同 地 理 分 
区 的 小 流域 ,以 NOAA-CPC-US 地 面 雨量 站 降水 再 
分 析 数 据 作为 参照 ,评估 不 同时 空 分 辨 率 的 卫星 
降水 产品 PERSIANN, PERSIANN- CDR, TRMM- 
3B42V7 , GPM-IMERG ,雷达 降水 StageIV 以 及 气候 
模式 降水 产品 ERA5 的 精度 ,并 结合 CREST 分 布 式 
水 文 模型 ,模拟 各 小 流域 的 日 径流 量 过 程 ,评估 不 
同时 空 分 辩 率 的 多 源 降 水 产品 在 小 流域 尺度 的 水 


文 模拟 效用 ,从 而 为 上 述 多 源 降 水 产品 在 水 文 行业 
中 研究 和 应 用 提供 参考 。 


2 研究 区 域 .数据 和 方法 


2.1 研究 区 域 

内 外 对 小 流域 的 概念 尚 没有 统一 的 定义 。 
不 同 国家 和 组 织 ,出 于 研究 目的 和 角度 的 不 同 , 划 
定 小 流域 的 标准 也 不 同 。 本 文 所 讨论 的 小 流域 ,是 
依据 美国 能 源 部 橡树 岭 国家 实验 室 (Oak Ridge Na- 
tional Laboratory, USA) 的 河流 分 级 系统 (USSCE， 
The U.S. Stream Classification System) 。 该 系统 将 美 
国 近 260 万 条 河流 分 为 8 级 ,从 大 到 小 分 别 为 超大 
河流 Great River, AY] Large River, 干流 Main stem, 
中 等 河流 Medium River, /]\ 7] Small River, 1 Vit 
Large Creek .小 溪 Creek ,河源 Headwater ? o ZI XC TE 
美国 大 陆 区 域 的 不 同 地 理 分 区 ,选取 了 9 个 闭合 小 
流域 。 这 些小 流域 中 的 河流 在 该 河流 分 级 系统 中 ， 
属于 河源 Headwater 或 小 溪 Creek 级别。 这 9 个 小 流 
域 的 地 理 分 区 ,流域 面积 ,年 均 流 量 等 信息 见 表 1。 
2.2 研究 数据 

本 文选 用 了 NOAA-CPC-US, PERSIANN ,PER- 
SIANN-CDR ,TRMM-3B42V7, GPM-IMERG, Stage- 
IV 和 ERA5 共 7 种 降水 数据 。 数 据 时 间 段 为 2002 年 
1 月 1 日 ~2018 年 6 月 30 日 。 研 究 中 设置 时 间 分 辨 
率 为 日 尺度 ,空间 分 辩 率 为 0.0083°*x 0.0083°( 约 
1km)。 由 于 上 述 7 种 降水 数据 的 时 空 分 辩 率 不 同 ， 
此 在 数据 预 处 理 过 程 中 ,将 时 间 分 辨 紊 小 于 日 尺 
度 的 数据 ,通过 累积 转换 为 日 尺度 ,空间 分 辨 率 则 
通过 双 线 性 插值 统一 到 0.0083° x 0.0083°。 
2.3 方法 
2.3.1 NOAA-CPC-US 降 水 


该 数据 集 是 美国 国家 


表 1 选取 的 9 个 美国 不 同 地 理 区 域 小 流域 的 基本 信息 


Tab.1 Nine small watersheds in different geographic regions of the United States 


序号 流域 编码 流域 名 称 河流 等 级 面积 /km? 地 理 分 区 年 均 流量 /mSec” 
1 HUCO0107000107 Smith River Headwater 221.51 Northeast(NE) 4.18 
2 HUCO0312000306 Upper Telogia Creek Headwater 348.89 Southeast(SE) 5.91 
3 HUC0602000211 Oostanaula Creek Headwater 182.00 Central 2.60 
4 HUCO0703000511 Kinnickinnic River Headwater 444.00 East North Central(ENC) 3.16 
5 HUC1012010909 Upper Spring Creek Headwater 519.00 West North Central(WNC) 0.61 
6 HUC1 113030208 Little washita River Headwater 625.85 South 1.48 
7 HUC1506010304 Cherry Creek Headwater 720.00 Southwest(SW) 0.91 
8 HUC1710010104 Calawah River Creek 351.88 Northwest(NW) 28.04 
9 HUC1801021102 Trinity River Creek 410.33 West 5.16 
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海洋 和 大 气 局 气候 预测 中 心 (NOAA Climate Predic- 
tion Center) 基 于 地 面 观测 降水 的 一 个 再 分 析 降 水 产 
品 。 该 数据 集 汇集 上 万 个 地 面 观测 台 站 降水 观测 
信息 ,并 与 雷达 ,卫星 观测 以 及 数值 模型 预测 的 降 
水 进行 比较 ,实施 了 质量 控制 。 在 考虑 地 形 效 应 
后 ,对 数据 进行 优化 插值 ,创建 日 降水 场 ,最 终 形成 
具有 更 高 定量 准确 性 和 一 致 性 的 降水 产品 。 其 空 
间 分 辨 率 为 0.25° , IE [E] Pel sic? E] 1970 4E 1 H 1 H 8 
今 ,空间 范围 为 20.0~50.0°N, 230.0 ~ 305.0 E, 4 it T 
美国 大 陆地 区 。 数 据 下 载 地 址 https://www.esrl.noaa. 
gov/psd/data/gridded/data.unified.html o 
2.3. NCEP-StageIV 雷达 降水 NCEP-StagelV 是 
美国 气象 局 国家 环境 预报 中 心 的 雷达 定量 估 测 降 
水 产品 ,项 目 起 始 于 2001 年 12 月 1 日 ,数据 分 发 起 
台 于 2002 年 1 月 1 日 。 其 数据 结合 全 美 大 约 150 个 
多 普 勒 新 一 代 天 气 雷 达 的 降水 估计 值 和 约 5 500 个 
地 面 雨量 计 的 测量 值 ,由 分 布 在 美国 大 陆 的 12 个 河 
流 预 报 中 心 按 各 自 区 域 生 产 。 数 据 采 集 开 始 于 每 
小 时 的 33 分 钟 ,并 以 1h 和 6h 为 时 间 步 长 进行 累积 
计算 ,迭代 复核 和 人 工 质量 控制 ,最 后 将 各 区 域 数 
据 馈 区 合成 美国 大 陆 雷 达 降 水 定量 估 测 数据 。 
该 数据 空间 分 辨 率 4 km, 时 间 分 辨 率 1h,6n,24h， 
降水 单位 mm。 数 据 详细 特性 及 下 载 链接 可 参考 
National Stage IV QPE Product 网 站 (www.emc.ncep. 
noaa.gov/mmb/ylin/pepanl/stage4/) 。 本 文选 用 Stage- 
TV 时 间 分 辩 率 6h 降 水 估 测 数据 。 
2.3.3 PERSIANN 降水 PERSIANN 降水 数据 集 由 
美国 加 州 大 学 欧文 分 校 (University of California, Ir- 
vine ) 水 文 气象 与 遥感 研究 中 心 (The Center for Hy- 
drometeorology and Remote Sensing, CHRS) 开 发 。 它 
基于 地 球 同 步 卫星 (GOES-8, GOES- 10, GMS-5, 
Metsat-6, Metsat-7) 提 供 的 红外 亮 温 图 像 ,使 用 人 工 
神经 网 络 估算 30 分 钟 降雨 率 , 然 后 汇总 到 1,3 ,6 小 
时 和 日 尺度 的 累积 降雨 量 , 并 使 用 低 轨道 卫星 
(TRMM, NOAA-15, -16, - 17, DMSP F13, F14, F15) 
降水 估算 产品 定期 修正 。 该 数据 产品 时 间 覆 羔 
2000 年 3 月 至 今 ,空间 范围 60°S ~ 60°N, 空 间 分 辨 
率 0.25°。 数 据 下 载 源 https://chrsdata.eng.uci.edu/。 
2.3.4 PERSIANN-CDR 降 水 ”该 数据 集 提供 了 从 
1983 年 至 今 ,纬度 范围 为 60°S ~ 60°N, 空 间 分 状 率 
23 0.25 ° ,时间 尺 度 为 日 的 降雨 量 估 算 值 。 它 由 美 
国 加 州 大 学 欧文 分 校 开 发 的 一 种 基于 卫星 还 感 数 


据 的 降水 产品 。 它 在 GridSat-B1 红 外 卫星 数据 上 通 
过 人 工 神 经 网 络 算法 (PERSIANN) 产 生 ; 人 工 神经 
网 络 算法 的 训练 是 用 美国 环境 预报 中 心 (NCEP) 
StagelV 雷达 降水 数据 完成 ;然后 ,使 用 全 球 降水 气 
候 学 项 目 (GPCP) 月 降水 产品 (GPCP v2.2) 进 行 调 
Re ,最 终 形成 该 数据 产品 。 美 家 研究 委员 会 
(NRC) 将 该 数据 定义 为 一 份 气候 数据 记录 (Climate 
Data Record) ,拥有 足够 长 度 ,一 致 性 和 连续 性 的 时 
间 序 列 ,用 以 确定 气候 变 率 和 变化 。 数 据 下 载 源 与 
PERSIANN 相同 。 

2.35 TRMM-3B42 V7 卫星 降水 TRMM 卫星 
2015 年 4 月 由 于 燃料 耗 尽 ,正式 结束 观测 任务 ,并 于 
同年 6 月 重新 进入 地 球 大 气 层 , 坠 人 南 印度 洋 。 由 
于 TMPA 降水 产品 的 成 功 ,TRMM 团队 并 没有 因为 
卫星 终止 观测 而 马上 停止 TMPA 降水 数据 产品 。 而 
是 ,运用 多 种 方法 来 延续 TMPA 降水 产品 ,计划 一 直 
更 新 到 2019 年 中 期 站 。 本 文选 择 TRMM-3B42 v7 
降水 产品 也 意 在 评估 当前 条 件 下 ,该 数据 集 模拟 水 
文 效用 的 能 力 。TRMM-3B42 V7 数据 空间 分 辩 率 
290.25 * ,时间 分 辩 率 为 3h ,降水 单位 为 mmh, Z 
数据 通过 NASA 的 TRMM 卫星 数据 网 站 下 载 (https:// 
pmm.nasa.gov/trmmy/)。 

2.3.6 GPM-IMERG 卫星 降水 ”GPM 是 美国 国家 
航空 航天 局 (National Aeronautics and Space Adminis- 
tration, NASA) 主 导 的 TRMM 的 后 续 全 球 卫 星 降 水 观 
测 计划 。GPM 的 核心 观测 平台 于 2014 年 2 月 28 日 
发 射 ,能够 提供 全 球 3h 以 内 雨 雪 观测 数据 。 相 较 
TRMM,CGPM 覆 盖 范 围 大 ,延伸 至 地 球 两 极 ;时 空 分 
辩 率 高 ,传感器 性 能 更 优 , 可 精确 的 探测 微量 降 
水 。 有 关 该 降水 产品 的 其 它 特 性 可 参考 其 门户 网 
站 提供 的 技术 文档 (https://pmm.gsfc.nasa.gov/GPM)， 
GPM 数据 的 Level 1 ~ 3 产品 也 可 通过 该 网 站 下 
载 。 本 文选 用 GPM 多 卫星 反 演 全 球 范围 30 分 钟 十 
雪 数 据 产品 (Integrated Multi-satellitE Retrievals for 
GPM,IMERG), 版 本 05B_Final。 该 降水 数据 通过 最 
新 GPROF2017 算 法 制作 ,校准 并 融合 了 了 GPM 星座 
全 部 微波 红外 以 及 地 面 观 测 和 其 他 可 能 的 传感器 
数据 Yi。 其 空间 分 状 率 为 0.1“, 时 间 分 辨 率 为 30 分 
钟 , 降 水 单位 为 mm:h ~o 

2.3.7 ERAS 再 分 析 数 据 集 降水 ”ERA5 是 欧洲 中 
期 天 气 预报 中 心 (ECMWF) 发 布 的 第 五 代 再 分 析 数 
据 集 。 相 较 于 ERA-Interim 和 ERA-40 前 代 再 分 析 
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资料 ,ERA5 采 用 了 改进 的 数据 同化 系统 (IFS Cycle 
41r2) ,四 维 变 分 同化 分 析 (4D-Var) 以 及 辐射 变 分 
偏差 校正 ,系统 性 地 提升 了 该 数据 集 的 质量 。ERA5 
提供 了 全 球 每 天 每 小 时 的 估计 值 ,而 且 每 小 时 输出 
分 辨 率 相对 于 ERA-Interim 而 言 有 相当 大 的 改进 。 
该 数据 包括 蒸发 ,湿度 ,气温 ,水 汽 压 ,降水 等 265 种 
变量 ,时 间 履 盖 范 围 1979 年 至 今 ,空间 分 辨 率 0.25” 
x 0.25" ,空间 覆盖 范围 89.142*S ~ 89.142°N, 180°W ~ 
180°E 。 本 文选 用 该 数据 集中 的 小 时 尺度 降水 数 
据 , 数 据 下 载 源 https://eds.climate.copernicus.eu。 
2.3.8 流域 径流 、 地 形 及 潜在 蒜 散 发 数据 ASCO 
个 小 流域 的 流域 边界 ,数字 地 形 DEM 和 径流 资料 ， 
分 别 采用 美国 地 质 调查 局 (USGS) 提 供 的 流域 边界 数 
Hi 4°! (Watershed Boundary Dataset V2.2.1) , 美国 
10m 精度 数字 高 程 模 型 (DEM) 以 及 美国 国家 水 信息 
系统 数据 。 在 DEM 数据 的 基础 上 ,应 用 GIS 空 间 分 
析 工 具 , 提 取 流 域 河 网 , 河 网 流向 , 河 网 汇聚 等 信 
息 。 本 文 用 分 布 式 水 文 模型 ,评估 各 降水 数据 水 文 
模拟 效用 时 ,需要 的 潜在 蒸 散发 数据 (PET) ,来 自 
USGS 饥荒 预警 系统 (Famine Early Warning System, 
FEWS)。 其 每 日 全 球 潜在 蒸 散 量 是 根据 每 6 小 时 从 
全 球 资料 同化 系统 提取 的 气候 数据 , 使 用 Penman- 
Monteith 方 程 计算 得 出 ,然后 求 和 获得 每 日 总 计 。 
该 数据 空间 分 辨 率 1.0° x 1.0 ,时 间 分 辨 率 1d。 数 
据 下 载 源 https://earlywarning.usgs.gov/fews/datadown- 
loads/Global/PET ; 
24 CREST 分布 式 水 文 模型 

CREST (Coupled Routing and Excess Storage ) 模 
73 rh KEREM Sy AA IK XC UAR 51 EL s 
(HyDROS) 和 美国 宇航 局 (NASA) SERVIR Jit H 
队 联 合 开 发 。 它 是 一 个 分 布 式 水 文 模型 ,将 全 流域 
划分 为 在 干 网 格 ,使 用 可 变 渗透 能 力 曲线 逐 网 格 进 
行 产 流 计 算 ,并 利用 多 线性 水 库 逐 网 格 计算 地 表 和 
地 下 水 汇流 ,最 后 通过 耦合 产 流 要 素 和 汇流 结构 , 
模拟 再 现 径流 过 程 ” 。 该 模型 在 全 球 尺 度 , 大 区 域 
尺度 ,以 及 中 小 尺度 都 有 出 色 的 水 文 模拟 和 预报 能 
力 ”””。CREST 模 型 运行 时 ,需要 输入 的 数据 包括 
数据 流域 信息 (如 流域 边界 ,DEM ,流向 等 ) ,降水 ， 
潜在 蒸 散 发 ,流域 出 口 流 量 ,初始 条 件 以 及 模型 参 
数 。 模 型 输出 结果 包括 土壤 含水 量 ,土壤 湿度 ,地 
面 径 流 等 。 为 节省 篇 幅 , 有 关 模 型 数据 预 处 理 , 模 
型 参数 取 值 范围 等 详细 信息 请 参考 美国 俄 克 拉 答 


KF HyDROS 实验 室 CREST 模 型 网 址 (http://hy- 
dro.ou.edu/research/crest/ ) 
2.5 统计 评估 指标 

本 文选 择 了 相对 偏差 (Bias) , 均 方 根 误差 
(RMSE) ,平均 绝对 误差 (MAE) ,相关 系数 (Correla- 
tion Coefficient, CC) ,4 种 统计 验证 指标 用 于 定量 评 
估 多 源 降水 产品 的 精度 。 采 用 标准 差 (Standard De- 
viation, SD) ,纳什 效率 系数 (NSCE), 相 对 偏差 (Bias) 
和 相关 系数 (CC) 评 估 各 降水 产品 的 水 文 模拟 效用 。 


> Sim,- $ Obs, 
Bias = | —— — —— —— 100 (1) 
RMSE = (2) 
zd n ! B 
MAE - a 245m. Obs, (3) 
SD= (4) 


Y (Obs, - Obs)(Sim, - Sim) 
CC=—2! (5) 
(Xo - Obs)’ X (Sim, - Sim) 


i=l i=1 


n 


> (Obs, - Sim) 
NSCE = | - ii — ——— (6) 
Y (Obs; - Obs) 


式 中 ;Sim 代 表 多 源 降水 佑 测 数据 或 CREST 模 型 模 
拟 出 的 径流 数据 ,0bs 代表 参考 降水 数据 或 径流 站 
流量 观测 数据 ;n 为 模拟 值 或 观测 值 总 数 ,i 为 第 i 个 
模拟 值 或 观测 值 。 当 4 个 评估 指标 Bias=0% , 
RMSE=0, CC=1, NSCE=1 代表 最 好 的 评 佑 结果。 为 
了 细致 地 定量 评估 各 降水 产品 的 水 文 模拟 效用 , 根 
据 Moriasi 等 人 的 工作 号 ,将 NSCE 值 分 为 4 个 区 间 ， 
代表 4 个 级 别 水 文 模拟 效用 : 较 差 (NSCE<0.5) , 适 
用 (0.50<NSCE<0.65) ,良好 (0.65 < NSCE<0.75) , 
优秀 (0.75 < NSCE<1.00)。 


3 结果 与 分 析 


3.1 降水 精度 评估 
由 于 降水 产品 的 降水 估算 精度 对 其 水 文 模拟 
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效用 有 较 大 影响 。 因 此 ,首先 对 各 降水 产品 的 降水 
精度 做 出 评估 。 因 为 NOAA-CPC-US 降水 数据 集 
汇集 了 上 万 个 地 面 观测 台 站 降水 观测 信息 ,然后 在 


考虑 地 形 效 应 后 ,对 数据 进行 优化 插值 ,并 实施 质 
量 控制 。 所 以 ,本文 以 2014 年 6 月 1 日 ~2018 年 6 


月 30 日 的 NOAA-CPC-US 降 水 产品 作为 参照 ,对 同 
时 期 的 其 它 降水 产品 StageIV , PERSIANN , PERSI- 
ANN-CDR , TRMM-3B42V7, GPM-IMERG , ERAS if 
行 降水 精度 评估 。 

3.1.1 要 格 尺 度 降水 精度 ”在 各 个 小 流域 ,以 1 km 
栅 格 尺度 ,计算 得 到 待 评估 降水 产品 同 NOAA- 
CPC-US 降水 的 相关 系数 (图 1, 图 2)。 在 6 种 降水 
产品 中 ,StagelV 与 NOAA-CPC-US 数 据 的 相关 系数 
最 高 ,平均 值 达 到 0.73 ,相关 系数 在 不 同 地 区 的 小 流 


域 的 变化 为 0.65 ~ 0.81; 在 West 和 Southwest 地 区 相 
关系 数 较 低 。 其 次 是 PERSIANN-CDR 和 GPM- 
IMERG 产品 ,与 NOAA-CPC-US 数 据 的 相关 系数 分 
别 为 0.51~ 0.69,0.53 ~ 0.73。 这 2 种 降水 产品 均 在 
Northwest 和 West 地 区 相关 系数 较 小 。 再 次 是 PER- 
SIANN 和 ERAS j* i , Ej NOAA-CPC-US 数据 的 相 
关系 数 分 别 为 0.49 ~ 0.66,0.49 ~ 0.62。 其 中 PERSL 
ANN 产 品 在 East North Central, West North Central 以 
及 Northwest 地 区 的 相关 系数 较 低 。ERA5 属于 气候 
模式 的 降水 产品 ,在 9 个 小 流域 的 平均 相关 系数 为 
0.55, 相关 系数 最 低 的 地 区 是 West, Northeast 和 
West North Central 地 区。 与 NOAA-CPC-US 数 据 的 
相关 系数 最 低 的 降水 产品 是 TRMM-3B42 V7, 相 关 
系数 在 0.44 ~ 0.59 ,在 9 个 小 流域 的 平均 相关 系数 只 


1.0 
0.8 上 I wl 
六 0.6 / 
dn 1012010909 E TI 
= 107000107 == 
= 312000306 = 1113030208 . 
| = 602000211 == 1506010304 
== 703000511 == 1710010104 ”最 人 人 
= 1801021102 
0.2 1 | fi 


StagelV PERSIANN 


| 
PERSIANN-CDR 


| | 
TRMM GPM ERAS 


1 多 源 降 水 产品 在 各 小 流域 与 NOAA-CPC-US 降水 相关 系数 


Fig. 1 Correlation coefficients between multi-source precipitation products and NOAA-CPC-US precipitation in each small watershed 


1.0 
s | ||! il | | 
ag 0.6 H | | | i 
l | I | | 
0.4 4 
um StagelV = TRMM - Max 
ums PERSIANN um GPM * Mean 
02 m= PERSIANN-CDR ma ERAS - Min 
Central ENC NE NW SE SW South WNC West 


图 2 不 同 地 区 多 源 降水 产品 与 NOAA-CPC-US 降水 相关 系数 


Fig. 2 Correlation coefficients between multi-source precipitation products in different regions and NOAA-CPC-US precipitation 
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^H 0.52, JE H YE Northeast, West North Central, East 
North Central FI Northwest 地 区 相关 系数 最 低 。 以 上 
6 种 降水 产品 与 NOAA-CPC-US 降水 的 相关 系数 均 
通过 0.05 置信 水 平 检验 。 

从 各 降水 产品 在 各 人 研究 区 域 相对 于 与 NOAA- 
CPC-US 降 水 的 均 方 根 误差 和 绝对 误差 ( 表 2) ,也 反 
映 出 降水 精度 的 差异 。 各 降水 产品 在 9 个 小 流域 的 
均 方 根 误差 ,绝对 误差 的 平均 值 ,StagelV 降水 RMSE 
为 5.78 mm, MAE 73 2.15 mm;PERSIANN-CDR 降水 
的 RMSE 为 6.52 mm, MAE 为 6.39 mm ; PERSIANN 
降水 的 RMSE 为 12.14 mm, MAE 6.90 mm; GPM- 
IMERG 卫星 降水 的 RMSE 为 23.27 mm,MAE 为 6.44 
mm; ERAS 降水 的 RMSE 为 53.02 mm, MAE 为 6.84 
mm; TRMM- 3B42 V7 卫星 降水 的 RMSE 为 65.92 
mm,MAE 为 7.06 mm。 

综 上 统计 指标 分 析 , 表 明 6 种 降水 产品 与 NO- 
AA-CPC-US 降 水 的 相关 性 有 差异 ,这 种 差异 具有 
地 区 特点 ,显示 出 降水 产品 的 空间 变异 性 。 降 水 产 
品 在 美国 高 纬度 ( 偏 北 ) 和 西部 山地 等 区 域 估 算 精 
度 略 低 。 各 降水 产品 与 NOAA-CPC-US 降水 吻合 
程度 从 高 到 低 ,依次 是 StageIV 雷达 降水 , PERSI- 
ANN-CDR fill GPM-IMERG 次 之 ,再 次 是 PERSIANN 
fil ERAS ,最 后 是 TRMM-3B42 V7。 

3.1.2 不 同 降水 强度 级 别 精 度 “ 本 文采 用 国际 气象 
组 织 (WMO ) 降 水 强度 级 别 分 类 标准 ” ,在 不 同 降水 
强度 等 级 上 ,更 精细 地 分 析 了 各 降水 产品 的 降水 精 
度 。 降 水 强度 按 每 小 时 降水 量 划 分 为 6 个 等 级 ,分 
别 是 (1) 无 十 :降水 <l mm; (2) 微 雨 :1 mm 降水 < 
2 mm; (3) 小 雨 :2 mm 降水 «5 mm; (4) PRY: 5 


mm 和 降水 <10 mm, (5) KEN: 10 mm ME 7K «20 
mm, (6) 28h :20 mm 二 降水 < 50 mm, (7) 大 暴雨 : 降 
水 =50 mm。 统 计 分 析 各 降水 强度 等 级 的 发 生 频 率 
(图 3) 和 各 降水 产品 在 不 同 降水 强度 级 别 的 频率 误 
差 ( 图 4)。 

相对 NOAA- CPC- US 降水 ,StageIV, PERSI- 
ANN, PERSIANN- CDR, TRMM- 3B42V7, GPM- 
IMERG, ERAS 降水 产品 在 不 同 地 区 9 个 小 流域 不 同 
降水 强度 的 发 生 频 率 并 不 相同 ,有 不 同 程度 上 的 高 
估 和 低估 ,误差 均 未 超过 +10 %。 总 体 上 ,多 数 降 水 
产品 在 无 十 和 微 雨 级 别 存在 低估 ,对 大 雨 和 暴雨 存 
在 高 估 。StageIV 雷达 降水 在 各 降水 强度 等 级 的 估 
测 精 度 都 高 于 PERSIANN，PERSIANN- CDR, 
TRMM-3B42V7 ,GPM-IMERG ,ERA5 降水 。 
3.2 水 文 模拟 效用 评估 

对 降水 精度 评估 之 后 ,本 文 用 上 述 7 种 降水 产 
品 分 别 驱 动 分 布 式 水 文 模 型 CREST ,模拟 美国 不 同 
地 区 9 个 小 流域 的 径流 过 程 ,评估 各 降水 产品 的 水 
文 模拟 效用 。 径 流 模拟 情景 设 定 2002 年 1 月 1 日 ~ 
2003 年 12 月 31 日 为 模型 预 热 期 ,2004 年 1 月 1 日 ~ 
2014 年 6 月 30 日 为 模型 率 定 期 ,2014 年 7 月 1 日 ~ 
2018 年 6 月 30 日 为 模型 验证 期 。 在 相同 的 率 定 期 
内 ,分 别 使 用 7 种 降水 产品 率 定 CREST 分布 式 水 文 
模型 参数 ,基于 各 自 降 水 产品 率 定 的 参数 集 ,然后 
在 相同 的 验证 期 对 径流 进行 模拟 。 采 用 NSCE, Bi- 
as, SD 和 CC 统计 指标 ,评估 水 文 模拟 效用 。 这 种 情 
景 设置 有 利于 独立 评估 各 降水 数据 的 流域 水 文 过 
程 模拟 效果 。 

由 图 5, 图 6, 图 7 分 别 说 明 各 降水 产品 在 美国 


表 2 多 源 降水 产品 在 各 研究 区 域 相对 于 与 NOAA-CPC-US 降水 的 误差 /mm 


Tab.2 Errors of multi-source precipitation products relative to NOAA-CPC-US precipitation in each small watershed / mm 


StagelV PERSIANN 


PERSIANN-CDR 


流域 


RMSE MAE RMSE MAE RMSE 


HUCO0107000107 5.63 1.89 10.22 — 6.95 5.82 
HUC0312000306 5.02 1.62 9.56 6.59 5.51 
HUC0602000211 4.97 1.66 1002 . 6.55 5.43 
HUCO0703000511 5.29 2.26 1823 7.46 6.27 
HUCIO12010909 4.99 2.03 19.68 7.61 6.54 
HUC1113030208 5.96 1.86 7.36 6.16 6.19 
HUC1506010304 6.58 2.34 8.36 6.14 7.29 
HUCI710010104 6.91 321 9.12 7.04 7.59 
HUC1801021102 — 6.71 2.46 1673 6.98 8.03 
均值 5.78 2.15 12.14 6.83 6.52 


TRMM-3B42V7 GPM-IMERG ERAS 
MAE RMSE MAE RMSE MAE RMSE MAE 
6.36 69.06 7.56 27.22 6.42 58.37 1.32 
6.44 64.77 6.89 15.19 6.13 51.01 6.75 
6.41 65.11 6.64 22.01 6.52 47.11 6.59 
6.05 68.66 7.35 27.04 6.64 50.78 6.89 
6.37 72.61 7.49 26.85 6.49 60.26 7.06 
5.89 65.56 6.42 18.72 5.64 46.27 6.27 
6.02 58.26 6.51 20.23 6.47 49.06 6.83 
7.02 63.89 7.63 24.36 6.89 57.31 TL1 
6.96 65.39 7.01 27.83 6.74 57.02 7.03 
6.39 65.92 7.06 23.21 6.44 53.02 6.87 
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& 60 = TRMM x 50 = TRMM s = TRMM 
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图 3 多 源 降水 产品 在 不 同 降水 强度 的 发 生 频率 分 布 


Fig. 3 Frequency of occurrence of each precipitation product at different precipitation intensities 
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4 多 源 降 水 产品 相对 NOAA-CPC-US 降水 在 各 降水 强度 级 别 的 频率 误差 


Fig. 4 Frequency error of multi-source precipitation products relative to NOAA-CPC-US precipitation in each precipitation intensity 


202101.00027v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


FHEWY 


红色 表示 率 定期 
[0^ 黑色 表示 验证 期 
+ NOAA-CPC-US 


Standard deviation 
Standard deviation 
N 
e 


红色 表示 率 定期 
黑色 表示 验证 期 黑色 表示 验证 期 
+ NOAA-CPC-US + NOAA-CPC-US 
© Stage IV 


o 
io 
Standard deviation 


o& 1 
20 30 Obs 10 20 30 40 Obs5 10 15 20 
HUC0107000107 HUCO0312000306 HUC0602000211 
= 红色 表示 率 定 期 红色 表示 率 定期 
红色 表示 率 定期 isi T 
Rost P ROMA CPUS NOAA CICUS 


é * NOAA-CPC-US IV 
0 0.10.2 O, TANN, 
2210.29 3 ey, PERSIANN-CDR 
20 TF 904, 9. un 
E E 
g^ 8 
8 [| 5 
© d 
$^ E 
à s E 
Obs 5 10 15 20 " Obs 5 10 15 Obs 5 10 15 
HUCO0703000511 HUC1012010909 HUC1113030208 
- 红色 表示 率 定期 示 
ee 黑色 表示 验证 期 LUPO | 
$ NOAA-CPC-US + NOAA-CPC-US * NOAA-CPC-US 
0 Q x PERS 
s 0.1 0.293 CO : PERSIANN. CDR 5 PERSIANN-CDR 
4 GPM-IMERG GPM-IMERG 
; v ERAS 
8 8 
5 Mey: 号 
E 10 Ox E 
© © 
5 g 
E ES 3 
E 0.95 8 
n a 
Rom 0.99 AV 
y 1 su i X A 
Obs 5 10 15 Obs 50 100 40 
HUC1506010304 HUC1710010104 HUC1801021102 


图 5 多 源 降水 产品 在 美国 不 同 地 区 9 个 小 流域 水 文 模拟 的 SD,RMSE 和 CC 统计 值 (图 中 Obs 表 示 实 测 径 流 的 统计 值 ) 


Fig.5 SD, RMSE, and CC statistics of hydrological simulations of multi-source precipitation products in nine small watersheds in 


different regions of the United States (Obs in the figure represents the statistics of measured runoff) 


不 同 地 区 9 个 小 流域 日 径流 量 率 定 和 模拟 的 结 
包括 相对 于 实测 径流 量 的 标准 差 SD , 15 75 AR V 2S 
RMSE ,相关 系数 CC, 纳什 效率 系数 NSCE 以 及 相对 
偏差 Bias。 

各 降水 产品 在 率 定期 径流 模拟 性 能 好 于 验证 
期 ,验证 期 的 各 个 统计 指标 均 有 衰减 。7 种 降水 产 
品 在 美国 不 同 地 理 分 区 小 流域 的 径流 模拟 性 能 不 
一 , 因 篇 幅 有 限 , 此 处 仅 以 NSCE 指标 对 各 降水 产品 
在 美国 不 同 地 理 分 区 小 流域 水 文 模拟 的 效果 展开 


说 明 , 其 它 指标 请 参见 对 应 图 表 。 根 据 Moriasi 等 人 
的 工作 中 ,将 NSCE 值 分 为 4 个 区 间 , 代 表 4 个 级 别 
水 文 模拟 效用 : 较 差 (NSCE<0.5) ,适用 (0.50< 
NSCE<0.65) ,良好 (0.65 < NSCES0.75) ,优秀 (0.75 < 
NSCPE<I.00)。 整 理 7 种 降水 产品 在 9 个 小 流域 率 定 
期 和 验证 期 的 水 文 模拟 效用 级 别 如 表 3 所 示 。 
整体 上 ,NOAA-CPC-US 和 StagelV 降水 是 7 种 
降水 产品 中 水 文 模拟 效果 最 好 的 2 种 。 在 美国 不 同 
地 理 分 区 的 小 流域 保持 了 较 好 的 径流 模拟 特性 的 ， 
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图 6 多 源 降水 产品 在 美国 不 同 地 区 9 个 小 流域 水 文 模拟 的 纳什 效率 系数 
Fig.6 Nash efficiency coefficient of hydrological simulation of multi-source precipitation products in nine 


small watersheds in different regions of the United States 
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图 7 多 源 降水 产品 在 美国 不 同 地 区 9 个 小 流域 水 文 模拟 的 相对 偏差 


Fig.7 Relative deviations of multi-source precipitation products in hydrological simulations of nine small 


watersheds in different regions of the United States 
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表 3 多 源 降水 产品 在 各 小 流域 径流 模拟 效用 级 别 


Tab.3 Performance of multi-source precipitation products in the simulation of runoff in each small watershed 
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NSCE 级 别 多 数 达 到 了 良好 以 上 ,能 够 重 现 流域 行 
洪 期 间 的 日 径流 量 过 程 ;GPM-IMERG 卫星 降水 和 
PERSIANN-CDR 次 之 ,在 美国 中 低 维度 地 区 的 小 流 
域 表现 尚 可 ,在 北部 ,西部 地 区 模拟 效果 较 差 。 
PERSIANN, ERAS 降水 的 径流 模拟 能 力 不 及 NO- 
AA-CPC-US 和 StagelV 降水 ,对 洪峰 流量 的 模拟 效 
果 则 存在 明显 波动 ,不够 准确 ,在 美国 中 部 ,南部 模 
拟 性 能 较 好 ,但 在 北部 高 纬度 地 区 和 美国 西部 , 西 
南部 地 区 较 差 。TRMM-3B42V7 可 能 由 于 传感器 精 
度 和 分 辨 率 的 缘故 ,在 多 个 小 流域 的 径流 模拟 不 
理想 。 


4 结论 


v 


通过 遥感 技术 获得 降水 估 测 数据 ,进而 驱动 水 
文 模型 ,模拟 和 预测 不 同 尺度 的 水 文 过 程 , 是 流域 
水 文 水 资源 研究 的 主要 趋 势 。 这 类 数据 对 洪涝 , 飓 
风 , 滑 坡 等 自然 灾害 的 预测 预报 极为 重要 。 同 大 尺 
度 流域 相 比 ,小 流域 是 研究 小 微 尺 度 水 文 水 资源 、 
水 环境 的 理想 对 象 , 也 是 探 明 小 微 水 文 水 资源 系统 
演变 规律 ,寻找 水 资源 可 持续 利用 方法 的 重要 
手段 。 

(1) 近年 来 ,国内 外 学 者 基于 遥感 降水 产品 开 
展 了 诸多 全 球 尺 度 ,洲际 尺度 ,大 流域 尺度 的 与 水 
文 水 资源 相关 内 容 的 研究 。 但 是 ,遥感 降水 估 测 数 
据 在 小 流域 的 研究 和 应 用 仍 不 充分 。 因 此 ,本 研究 
以 NOAA-CPC- US 降水 作为 参照 ,对 比分 析 了 
StagelV, PERSIANN, PERSIANN- CDR, GPM- 
IMERG, ERAS, TRMM 3B42 V7 降水 的 精度 ,并 用 这 
7 种 降水 产品 驱动 CREST 分 布 式 水 文 模型 ,评估 了 


它们 在 小 流域 水 文 模拟 的 可 行 性 。 

(2) 从 降水 精度 评估 结果 来 看 ,各 降水 产品 在 
美国 北部 高 纬度 和 西部 山地 等 区 域 估算 精度 略 低 ， 
在 中 部 ,南部 ,东部 保持 较 好 的 精度 。 各 降水 产品 
与 NOAA-CPC-US 降 水 吻合 程度 从 高 到 低 , 依 次 是 
StagelV 雷达 降水 ,PERSIANN-CDR 和 GPM-IMERG 
次 之 ,再 次 是 PERSIANN 和 ERA5S, 最 后 是 TRMM- 
3B42V7。 从 对 不 同 降 水 强度 的 捕捉 能 力 来 看 , PER- 
SIANN, PERSIANN- CDR, GPM- IMERG 和 TRMM- 
3B42V7 卫星 降水 产品 在 无 南 、 微 雨 存 在 不 同 程度 的 
低估 ,大 雨 和 暴雨 存在 高 估 。 

(3) 在 水 文 模拟 效用 评 佑 中, 总体 上 NOAA- 
CPC-US Fil Stage IV 雷达 降水 在 各 小 流域 水 文 模拟 
中 效果 较 好 ,PERSIANN-CDR,CGPM-IMERG PER- 
SIANN, ERAS 降水 的 径流 模拟 能 力 次 之 ,TRMM- 
3B42V7 模 拟 效果 不 理想 。 基 于 本 文 设 置 的 水 文 模 
拟 情景 。 

上 述 结论 说 明 在 美国 NOAA-CPC-US 降水 和 
StagelV 雷达 降水 基本 能 满足 小 流域 水 文 模拟 的 需 
求 ;而 使 用 PERSIANN, PERSIANN- CDR, GPM- 
IMERG, ,降水 进行 小 流域 水 文 模拟 时 需要 审慎 , 尤 
其 在 在 美国 北部 ,西部 地 区 ;ERA5,TRMM-3B42V7 
降水 数据 尚 不 足以 支撑 小 流域 水 文 模 拟 。 
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Evaluating runoff simulation of multi-source precipitation data 
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Abstract: Small watersheds are ideal objects for studying the evolution of small and microscale hydrology and 
water resources systems. A small watershed is the smallest unit for calculating river water and sediment production 
and is the best regional scale for hydrology and soil erosion research and management. Through RS technology, a 
climate model obtains precipitation estimation data and drives a distributed hydrological model to simulate and 
predict hydrological processes, which clarifies the inevitable trends of basin hydrology and water resources 
research. Using NOAA- CPC- US precipitation data as a reference, an analysis of the PERSIANN, PERSIANN- 
CDR, TRMM-3B42V7, GPM- IMERG, StagelV, and ERAS precipitation data products were compared for nine 
small watersheds in different regions of the United States. The accuracy of theseseven precipitation products 
allowed them to drive the CREST distributed hydrological model, which evaluated the hydrological simulation 
effects of the precipitation products. The study shows that the NOAA- CPC- US precipitation data product is the 
highest, followed in decreasing order by StageIV, PERSIANN-CDR,GPM-IMERG, PERSIANN,ERA5, and TRMM- 
3B42V7. The precipitation estimation accuracy of each precipitation product in the small watersheds in the high 
latitudes and western mountains of the northern United States is lower; however, there is better precipitation 
accuracy in small watersheds in the central, southern, and eastern parts of the United States. In the hydrological 
simulation utility evaluation, the CREST model parameters were determined using seven kinds of precipitation 
products, respectively. After obtaining the set of parameters, the daily runoff process of the basin was simulated for 
the same verification period. The comparison results show that NOAA-CPC-US and StagelV have better effects on 
the hydrological simulation of small watersheds. However, caution should be exercised in hydrological simulations 
in the northern and western parts of the United States using PERSIANN, PERSIANN- CDR, GPM-IMERG, and 
ERAS precipitation data, and the TRMM—-3B42V7 simulation effect is not ideal. 
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